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APPLICATIONS OF ANALYTIC FUNCTION THEORY TO ANALYSIS 
OF SINGLE -SIDEBAND ANGLE -MODULATED SYSTEMS 

By John H. Painter 
Langley Research Center 

SUMMARY 

This paper applies the theory and notation of complex analytic time functions and 
stochastic processes to the investigation of the single -sideband angle -modulation process. 
Low -deviation modulation and linear product detection in the presence of noise are care- 
fully examined for the case of sinusoidal modulation. Modulation by an arbitrary number 
of sinusoids or by modulated subcarriers is considered. Also treated is a method for 
increasing modulation efficiency. The paper concludes with an examination of the asymp- 
totic (low-noise) performance of nonlinear detection of a single -sideband carrier which 
is heavily frequency modulated by a Gaussian process. 

INTRODUCTION 

The theory of complex analytic time functions and stochastic processes is applied 
herein to the mathematical examination of single - sideband (SSB) angle -modulated sinus- 
oidal carriers. Emphasis is on application of the theory, rather than on derivation of it. 
Elements of the theory and the accompanying notation are fully developed in reference 1. 

There are two purposes for this paper. First, and most important, is the desire 
to determine mathematically the communication efficiency of SSB angle -modulation. This 
determination is made separately for low-index sinusoidal modulation with product demod- 
ulation and for high-index random modulation with nonlinear (discriminator) demodula- 
tion. A secondary purpose is to investigate the application of the theory of analytic func- 
tions to the solution of SSB modulation and detection problems. 

First the basic mathematical models for band-pass noise and signals and for their 
correlation and spectral functions are established on the basis of the complex analytic 
notation. Next, "modulation functions" are defined for the signals to be treated, and low- 
deviation modulation by subcarriers is discussed. 

In the section on "Sinusoidal Single-Sideband Amplitude -Phase Modulation and 
Detection," four unique models for signals having one subcarrier are presented and their 
power structure is derived. An observation is made on the transmission efficiency 



achievable with the four signal models. Then a more efficient signal is synthesized from 
two of these models. One of the four models is examined and some observations are made 
for the case of an arbitrary number of modulating sinusoids. Then the consequences of 
modulated subcarriers are investigated. The section concludes with a modeling of linear 
product demodulation of one of the four signals in the presence of noise. 

The section on "Gaussian Single -Sideband Frequency Modulation" begins with 
derivations of the average power of a Gaussian SSB-FM signal. Next the equivalent 
phase of the sum of the modulated signal plus Gaussian channel noise is investigated. A 
low-noise approximation to the equivalent phase process is used to obtain an asymptotic 
result for the phase noise autocorrelation function. The section concludes with a 
modeling of the nonlinear demodulation of the signal in the presence of noise under the 
assumptions of low noise and high modulation. 

SYMBOLS AND NOTATION 

The symbols and mathematical notation for the body of the paper are defined below. 
Symbols used only in the appendixes are not listed here. Since the appendixes are short 
and self-contained, each new symbol is defined within each appendix. 

Subscripting, as applied to functions, is standard. Thus, the symbol ) means 

the function R( ) which is related to functions x( ) and y( ). In some cases double 
subscripting is used. Thus, the symbol P s ^ indicates the ith value of P g . General 
mathematical notation is standard except as detailed below. 

Mathematical notation: 


E{) 

statistical expectation 

j=\A 


L( ) 

linear transformation 

Re( ) 

real part 

( ) 

Hilbert transform 

1 1 

absolute value 

A 

equality by definition 
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approximate equality 


( )* complex conjugate 

( ) time derivative 

( )" second derivative 

Real numbers : 

A amplitude constant, dimensionless 

B bandwidth constant, hertz 

c arbitrary constant, real or complex, dimensionless 

N noise power constant, watts 

P, S signal component power, watts 

Q, ri noise spectral density amplitude, watts per radian per second 

T limit of integration, seconds 

t, t time variable, seconds 

j3 amplitude constant, radians 

0 constant or dummy variable, radians 

cr^ variance of a stochastic process, radians^ 

0 constant or dummy variable, radians 

co frequency variable or constant, radians per second 



Deterministic functions and stochastic processes: 


A( ) real envelope function of modulated signal (deterministic or stochastic, 

depending on context) 

f( ) real low-pass message function (deterministic or stochastic, depending on 

context) 

g( ) real low-pass message function (deterministic or stochastic, depending on 

context) 

I( ) modified Bessel function 

J( ) Bessel function 

m( ) complex low-pass modulation function (deterministic or stochastic, depending 

on context) 

N( ) noise spectral density 

n( ) real band-pass stochastic noise process 

R( ) correlation function or envelope function of signal plus noise, depending on 

context 

r( ) sum function of signal plus noise 

S( ) general spectral density (Fourier transform of R( )) 

s( ) real band-pass signal function (deterministic or stochastic, depending on 

context) 

x( ) real low-pass signal or noise function (deterministic or stochastic, depending 

on context) 

y( ) real low-pass signal or noise function (deterministic or stochastic, depending 

on context) 
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z( ) complex low -pass signal or noise function (deterministic or stochastic, 

depending on context) 

5( ) the delta function 

0( ) subcarrier phase function or noise perturbation of the equivalent phase func- 

tion of signal plus noise, depending on context 

A( ) first zone output function of an ideal limiter 

v( ) complex band-pass stochastic noise process 

<3?( ) statistical characteristic function 

4>( ) real phase function of a modulated signal (deterministic or stochastic, 

depending on context) 

\(/( ) complex band-pass signal function (deterministic or stochastic, depending 

on context) 

ANALYTIC SIGNAL AND NOISE NOTATION 

In the analysis of SSB angle -modulated systems, maximum use is made of analytic 
notation. That is, real signal and noise processes are expressed as the real parts of 
complex analytic processes to exploit the advantages inherent in such notation. Deriva- 
tions of the properties of the analytic notation have been made elsewhere (refs. 1 to 6). 
Here, the properties of the notation are simply assumed. 

Given a real time function x(t), which may be either a deterministic function 
(signal) or sample function from a random process (noise or noise -modulated carrier), 
x(t) may be represented as the real part of a complex analytic function z(t), where 

z(t) = x(t) + jy(t) (1) 

In the random case, it is assumed without loss of generality that z(t), x(t), and y(t) 
are zero -mean and weakly stationary. 

The modulated signal and noise processes dealt with in this paper are of the band- 
pass type. That is, the frequency spectra of both the signals and the accompanying noise 
are assumed to be band-limited. To describe such processes suitably it is convenient to 
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write them in "polar form" as z(t) exp(jco c t). The exponential function serves merely 
to translate the spectrum of z(t) to the neighborhood of the carrier frequency oo c . If 
z(t) is assumed to have a low -pass spectrum, then z(t) exp(jw c t) has a band-pass 
spectrum. Hence, in this paper it is assumed that z(t) is low -pass. Such a restric- 
tion is not necessary in general and is used strictly for simplification. 

Band-Pass Noise Processes 
A real band-pass noise process n(t) -is denoted as 

n (t) = Re[^(t)] (2) 

where 

v(t) = z(t) exp(ju> c t) (3) 


Thus, 


n(t) = x(t) cos w c t - y(t) sin co c t 


(4) 


Equation (4) is a well-known formulation for a band-pass noise process in terms of two 
low-pass processes, x(t) and y(t), and a reference frequency u> c (ref. 7). The 
requirements on x(t) and y(t) to insure that n(t) is zero-mean and weakly stationary 
are given in terms of the correlation functions as 


Rxx(t) — RyyC 7 ") 

RyxW = _ ^yx^ _T ) 


(5) 


where x(t) and y(t) are individually zero-mean. 

If cu c is in the set of frequencies for which the spectrum of n(t) is nonzero, 
then equations (5) are satisfied if x(t) and y(t) are individually and jointly Gaussian 
(ref. 8). If o> c is in the set of frequencies for which the spectrum of n(t) is zero, 
then x(t) + jy(t) is analytic, or conjugate analytic, and equations (5) are satisfied, as is 
the additional relation 


Ryx( T ) _ ±R-xx( T ) 


( 6 ) 
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where the caret denotes Hilbert transform (ref. 1). The plus and minus signs are chosen, 
respectively, according to whether the frequency domain of the spectrum of n(t) is above 
or below u> c . 

In dealing with signals which have single sidebands with respect to a carrier fre- 
quency w c , it is convenient to reference the channel noise to the same frequency o> c . 

In this case, the noise process is also single -sideband in nature and equation (6) applies. 
Although in this case it is not strictly necessary, the channel noise is taken as Gaussian. 

For computational use in the remainder of the paper, the relations between the 
spectral densities of x(t), y(t), and n(t) are now derived. An idealizing assumption is 

made for computational simplicity. It is assumed that the noise process n(t) results 
from filtering a white Gaussian process with an ideal upper sideband filter. That is, the 
spectral density Snnfw) of n(t) is taken as 


SnnM = 




(all other 


w) 


( 7 ) 


where c o c , the carrier reference frequency, is the "lower cutoff frequency" for S nn (u>) 
and <x> c + u> u is the "upper cutoff frequency." 

From equation (All) of appendix A, the spectra of x(t) and y(t) are 


SxxM = SyyM = 


fQ; (M = w u) 

0; (all other lo) 


( 8 ) 


Figure 1 is given to clarify the relations expressed in equations (7) and (8). 

s nnM 

i — i i 1 — Q 


-u>c 

0 

w c w c + 

S xx (a;) - Syy(cu) 



Q 

« 

* 



-Wq ”^U 8 W-jj U^c 

Figure 1.- Spectral relations. 
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Band-Pass Signal Processes 

A real band-pass signal process s(t) is handled in a manner similar to that for 
noise processes. That is, s(t) is taken as the real part of some analytic i//(t), where 
i//(t) has the same form as the v(t) in equation (3). That is, let 


s(t) = Re [^(t)] 
i^(t) = m(t). exp(jw c t) 


(9) 


where m(t) is a "modulation function" which is, in general, complex. 

The properties of m(t) may be determined from examination of s(t). The most 
general form for a "modulated carrier" is 


s(t) = A(t) cos|o> c t + <p( tjj (10) 

where A(t) is an amplitude function and <p(t) is a phase function. Various forms of 
modulation are obtained for different choices of A(t) and <p( t). For 0(t) held con- 
stant, double -sideband amplitude modulation results. For A(t) held constant, double- 
sideband phase or frequency modulation results. If A(t) and cp( t) vary simultaneously 
and if they are related functionally, then various hybrid modulations are obtainable which 
have, in general, nonsymmetric sidebands. 

To obtain equation (10) from equations (9) requires 


m(t) = A(t) exp[j<£(t)] (11) 

Now m(t) may be either deterministic or stochastic, depending on which type of modu- 
lated carrier is being considered. For the stochastic case, sufficient conditions for s(t) 

to be weakly stationary are for m(t) to be zero-mean and (from ref. 1): 

E {m(t+r) m(t)} = 0 (12) 


It should be noted in passing that for the stochastic case, the phase function <£( t) of m(t) 
may require a uniformly distributed random constant to satisfy equation (12). This matter 
is treated in detail in appendix B. 

Given the modulation function m(t), the autocorrelation and spectral density func- 
tions of the real signal s(t), in the stochastic case, are obtained from equations (A16) 
and (A22), respectively, as 
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r ss( t ) = | Re [ R mm( T ) exp(jco c T)J 
Sss( w ) 1 |'[pmm( a)_a, c) + Smm( a,+W c)| 


(13) 

(14) 


Modulation Functions 


Bedrosian (ref. 6) has given the modulation functions m(t) for many types of mod- 
ulated signals. Those of chief interest in this paper are the angle -modulated signals. By 
an angle -modulated signal, what is explicitly meant is a modulated carrier where a 
"message function" f(t) is linearly related to the phase function <p(t). That is, for 
some f(t) which is to be transmitted and recovered, f(t) must be recoverable by a 
linear operation on (p(t). There are also certain hybrid signals where f(t) is recov- 
ered linearly from s(t). 

For a deterministic message function f(t), the most general modulation func- 
tion m(t) which results in an SSB angle -modulated signal is 


m(t) = exp jc L 


[f(t)+jf(t) J 


(15) 


where L denotes a linear transformation and c is a complex constant. Here m(t) 
is the exponent function of an analytic function and, hence, is analytic. 

For stochastic s(t), it is sufficient for this analysis to consider m(t) as simply 


m(t) - exp 



exp(j e) 


(16) 


where f(t) + jg(t) is analytic, and 6 is uniformly distributed in n}. From appen- 

dix B it is seen that m(t), as given in equation (16), satisfies the stationarity require- 
ment of equation (12) regardless of the statistics of f(t) + jg(t). 

For the stationary m(t) of equation (16), the autocorrelation function is 


R mm( T ) = E (m(t+T) m*(t)} 

= E 

A *(l,j,-l,j;r) (17) 



+ jg(t+r) - f(t) + jg(t)1 


where <3> denotes characteristic function (ref. 9). 
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SINUSOIDAL SINGLE-SIDEBAND AMPLITUDE -PHASE 
MODULATION AND DETECTION 


Unique Modulated Carrier Forms 

From equation (15) may be written four unique modulation functions for SSB (upper - 
sideband) sinusoidal modulation: 


m x (t) = expjjjf(t) + jf(t)]J 
m 2 (t) = exp^-j|f(t) + jf(t)]| 

m 3 (t) = exp^l|f(t) + jf(t)]| 

m 4 (t) = exp^-ljf (t) + jf(t)]j 


(18) 


Now, denote 


f (t) ^ /3 cos w m t 


(19) 


where /3 and co m are constants. Then the analytic signals i//.(t), corresponding to 
the mi(t), for i = 1, 2, 3, and 4 are 


i^(t) = expj^/3 exp j(w m t + |j" j exp(jo) c t) 


'"m* + f 


exp(jo) c t) 


yield 


= exp |^-/3 expjj ^ 

!^ 3 (t) = exp^/3 exp|jw m t|| exp(jw c t) 

i^ 4 (t) = exp^-/3 expjj oj m tj| exp(jw c t) 

If the Maclaurin series expansion for the exponent function is used, equations (20) 


( 20 ) 
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J fr ex p(j [( w c + nw m) t+ n §f} 

n=0 ’ ^ 

^2^)= £ (- 1 ) n ^ ex p[j[( w c + na, m) t+n |]J 


c° ^ 

^3 (t)= I ^ exp { j [( 

r\ ^ 


w c + na >m 


°° n ^ 
^ 4 (t) = y (-l) n |r exp/j ( 


w c + nw m 


Now the real signals s^(t) corresponding to the ^(t) are given in the form of equa- 
tion (10) and in series form as 

°° > 

s^(t) = exp (-/3 sin w m t)cos (a> c t + /3 cos w m t) = ^ cos fru> c + nu> m )t + n ^1 

n=0 

00 n 

SgCt) = exp^/? sin w m t^cos ^w c t - (3 cos a> m t ) = S (- 1 ) n fr C0S [( w c + nw m)t+n|J ! 


Sg(t) = exp(/3cos w m t)cos (u> c t + /3 sin w m t) = cos 


a) c + na) m 


OO 

s 4 (t) = exp(-/3 cos a> m t)cos(w c t - /3 sin a> m t) = ^ (-l) n ^ cos 


a) c + nw m 


Although the autocorrelation functions of the physical signals are not of prime inter- 
est, the mean-squared values, or powers, are. In the deterministic case, the autocorre- 
lation function of the ith analytic signal may be defined as a time average: 


I Vi W ' 2? I.t 'M t+T) dt 
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The autocorrelation function of the ith physical signal is then 


Rc 


W T >4 Ite [ , w1] 

The power, or mean- squared value, of the ith physical signal is 


Thus 


P s . = £ Re 

s i 2 


P Si “ R SiS^°) 


T 

lim -L f ^(t) ^*(t) dt 

T—oo s i J .q 1 


Due to the symmetry properties of R^.^(t), equation (26) reduces to 

T 

P Si = I f ^i(t) ^i*(t) dt 

1 2 T— 00 2T A 

Now, it is easily verified that 

""exp (-2/3 sin w m t) ; (i = 1) 
exp(2/3 sin o> m t); (i = 2) 
exp(2/3 cos o) m t); (i = 3) 
exp (-2/3 cos w m t); (i = 4) 

Thus i/^(t) ^ (t) is periodic in 2ir for 0 m , a dummy variable defined as 


^i(t) = ( 




^m _ w m* 


It follows that 


?s i = hfio exp ( 2/3 C ° S 0m ) d0m = I (i = *’ 2; 3 ’ 4) 

where Iq is the modified Bessel function (ref. 10). 


( 24 ) 


(25) 


(26) 


(27) 


(28) 


(29) 


(30) 
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More Efficient Modulation 


From equations (22) it is seen that the power residing in the recoverable signal 
component, defined here as the first-order sideband, is 


P 



( 31 ) 


Thus, from equations (30) and (31) the ratio of recoverable sideband power to total power 
is 



U2/3) 


(32) 


This function of modulation index /3 has an absolute maximum, 



0.475; (/3 = 1.29) 


(33) 


If efficiency is defined as percentage of recoverable signal power, referenced to total 
power, then this modulation process has a maximum efficiency of 47.5 percent. It is 
desirable to ascertain if there exist methods of increasing this efficiency. 

Figure 2 shows in three dimensions the phase relationship for the residual carrier 
and sideband components for the four signals of equations (22). It is easily determined 
from equations (22) and (30) that for the most efficient case, the sum of the powers in the 
residual carrier and the first- and second-order sidebands accounts for 99 percent of the 
total power. Thus, if two of the four signals can be combined so that the first-order 
components add and the carrier and second-order components subtract, then the overall 
efficiency can be increased. Figure 2 shows that this can be accomplished by subtracting 
S2(t) from Sj(t) or s^(t) from Sg(t). Furthermore, by adjusting the amplitudes of 
the two signals, the relative amplitudes of the carrier and first-order sideband can be 
controlled. Finally, by adjusting the modulation indices, the relative amplitudes of the 
carrier and second-order sideband can be controlled. 


Although it is not the purpose of this paper to give a complete synthesis technique 
(ref. 11) for sinusoidally modulated SSB-AM-FM signals, it is interesting to see, for the 
signal constructed above, how solutions for j3 can be formulated. Define the combina- 
tion signal as 


s(t) = s 2 (t) - s x (t) 


(34) 
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+ CcJj-q Ct ) q + 2 CJj-Q Cl)q 4- 3 Cl) 2 


/ 



Figure 2.- Relative phasing. 

where 


s-^t) = exp^-/3j sin co m t) cos(cu c t + jSj cos o> m t) 
s 2 (t) = A 2 exp(j3 2 sin c<; m t) cos(o> c t - /3 2 cos w m t) 


From equations (22) it can be determined that the power of the residual carrier is 


p c = \ ( A 2 - A l) : 


Likewise the powers of the first- and second-order sidebands are 


P 1 - |( A 2^2 + A l ^ 1 ) 2 

V^M-Vi 2 J 
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Now define 


Suppose constants k^ 


Then 


Equations (40) may be solved implicitly for C a and in terms of design constants 
kj and k2 with /3j as parameter. Then may be varied to satisfy some other 
constraint, say, on the total power P s . 

An expression for the total power P s in s(t) may be obtained in the same 
way P Si was obtained for the single signal. Let the analytic signal i^ s (t) be 

(41) 


(42) 


(43) 


where 


and 


m, 


m 


^ s (t) = ^2 (t) " ^j(t) 

i// 2 (t) = A 2 m 2 (t) exp(joj c t) 
^j(t) = Ajm^t) exp(jco c t) 

>(t) = exp |jj[f 2 (t) + jf 2 (t)J J 
L (t) = exp|jjf j(t) + jf^t)] j 


^ A Q a C a 

A! = C “ h - C t> 

and k 2 are defined as 


( 38 ) 


A k 2 

Pc = 1 


2 

p = k 2 


(39) 


k JjSE. <£lLH 


Cor-l 


k 2 ~ 


tfjcoFis 2 - x ) 

2 (C a - 1) 


) (o^l) 


(40) 
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where 



fi(t) = ^ cos w m t 1 

fi(t) = & sin o) m t J 

(44) 

Then 


p s (t) = jA 2 m 2 (t) - A 1 m 1 (tj] exp(jo; c t) 

(45) 

and 

P s (t) P s *(t) 

.= |A 2 m 2 (t) - Ajm^tj] [l 2 m 2 *(t) - A 1 m 1 *(t)] 



= A 2 2 exp|2f 2 (t)j + Aj 2 exp |^2f^(t)j - 2A^A 2 expjf 2 (t) - fj(t)j cosjf 
= A 2 exp(2/3 2 sin o> m tj + Aj exp(-2/3^ sin o> m tj 

2 (t) + f i (t)j 


- 2AjA 2 exp^ 2 - /3.|Jsin ^ m tjcos j(/3 2 + fi^jcos o> m tj 

(46) 


The total power is defined as a time average, as was done in equation (27): 

P s = I 2T It ^ * <47> 

Making the substitution as in equation (29), ^(^ni) ^ s *(0 m ) is periodic in the 

0 m interval - ^ . Thus 

.. •. /^37t/2 

P s=^J.,/2 * s ' *<e m ) <«> 

Now define the transformation 

* - «m - | < 49 > 

Then, advantage may be taken of the even symmetry of p s (4>) P s *(4>) in the Q-77, -tt] 
interval to define 
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By applying several Bessel function identities, P g is obtained as 



* A i%( 2 ' 3 i) - 


2A 1 A 2 


(/3 2 +/3i) I 0 (^2"^l) 


+ 


2 J J 2k( /3 2 + ^l) I 2k( /3 2“ /3 l) 
k=l 


(51) 


Substituting the definitions of equations (38) yields 


P s = t( C « 2i o( 2 7i) +I o( 2 M - 2C a { J o[( C /3 +1 )^l] I o[( C /3- 1 )%] 


OO 

+ 2 I 

k= 


, J 2k[( C g +1 )^ I 2k[( C /3- 1 )' i l]} 


(52) 


Thus the synthesis procedure for the combined SSB-PM signal is as follows: 

Solve equations (40) for C a and with parameter jSj, so that the kj and ^ 

constraints are satisfied. Vary the solution with /3^ so that the constraint on total 

power P_ is satisfied, 
s 


Modulation by Multiple Sinusoids 

So far, this section has dealt with modulation by a single sinusoid. The treatment 
is now generalized to an arbitrary number k of sinusoids. Let the modulation func- 
tion m(t) be 



m(t) = exp/ j 


y Pi exp^j^t) 


i=l 





(53) 


where ^ and are the modulation index and frequency, respectively, of the ith sub- 
carrier. Equation (53) shows that the complex modulation function for k sinusoids of 
frequency o^, for i = 1, 2, . . . k, may be expressed as the product of k infinite 
series. 

The Cauchy rule (ref. 12) for multiplication of two infinite series is 


OO 

“ 


OO 

l 

rv. 

m j 

X 

I la -2 

O 

II 

£ 



m 2 =° 

1— * 

— 




l 

2 -1 


I 

H= 




(54) 


where the superscript numbers simply serve to identify the terms of each series. This 
product rule may be extended to k-fold multiplication as 


k 

TTI 

i=l 


oo 

L V 


m^O 


n k 


= Z E 

n. — H n. ^ =l 


n k _n k-l 



n k n k-l"*** _n 2 


I 


^k^k-l 


k-1 



n, 


k-i 


( 55 ) 
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Substituting equation (53) into equation (55) gives 


m (t) = y ••• ^ 


n k~ n k-l •• • n 2/ . N n k" n k-l -, **“ n l 


n k =0 n k-r° 


x exp^j 


^0 \( n k - n k-l - - “ n l) 1 
k-1 

L 

i=l 


k- 1 (/3i) n k-i 


TTP 


iU ( n k-i)' 


K- JL 

(v + |)( n k - n k-l - - - n l) + JZ + |) n k-i 


( 56 ) 


Multiplying m(t) in equation (56) by exp(jw c t) and taking the real part of the 
product gives the resulting physical signal s(t): 


s(t) = 


TT ex P ("% sin "i 4 ) 

i=l 


cos w„t+ / 8, cos co.t) 

\ i=l / 


■L H ••• I 


n k -n k-r*”' n 2 


(ffk) 


n k n k-l •" n l 


n k =0 n k-l =0 


n^=0 


X COS 


( n k - n k-l - - - n i) : 
k-1 






“c‘ + ("t - n t-i - - - n i)(“k l + i) + H n k-i (“i‘ + 1) 

' ' i= 1 


(57) 


From equation (57), any arbitrary modulation term, say the jth, may be isolated. 
The jth term is explicitly the one for which the frequency is a> c + Wj and for which 



Denoting the jth modulation term as Sj(t), equations (57) and (58) yield 


Sj (t) = /3j cos [(Wc+ w j) t+ f 


(58) 


(59) 
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Equation (59) shows the major difference between the power behavior of this hybrid 
modulation and classical angle modulation. In the usual double -sideband angle -modulated 
signal, the amplitude of the jth modulation term is a product of Bessel functions whose 
arguments include all the "modulation indices” ^ (ref. 13). Thus, the amplitude of the 
jth term is reduced by the presence of the other modulating sinusoids so that the total 
signal power is constant and is just equal to the unmodulated carrier power. In this 
hybrid SSB modulation scheme, the amplitude of the jth term is a function of /3j only. 
Hence, the total signal power increases with the addition of modulating sinusoids. This 
behavior is attributable to the amplitude modulation present in the hybrid process. 


Subcarrier Modulation 

As in classic double -sideband modulation schemes, it is possible to modulate the 
modulating sinusoids to obtain subcarriers which then are processed as in the previous 
section to obtain an SSB spectrum. In general the subcarriers may be amplitude modu- 
lated, angle modulated, or hybrid modulated to obtain subcarriers which, themselves, 
have SSB spectra. However, the following development is limited to classically angle- 
modulated subcarriers. 

Now, modify the message function f(t) of equation (19) to be 


f(t) = j3 cosjw m t + 0(t)j 


(60) 


where 0(t) represents the phase function of the subcarrier. It is assumed that 0(t) 
is linearly related to some "submessage function" of f(t). The Hilbert transform 
of f(t) is given as 


f(t) = /3 sin [co m t + 0(tj] 


(61) 


provided the spectrum of exp|j0(tj| is zero for co < -o> m (ref. 14). Then it is easily 
seen that the four signals of equations (22) are modified to 

Sj(t) = exp ^-/3 sin [w m t + 0(t)jj cos ju> c t + /3 cos[o> m t + 0(t)jj 


oo 

= ^ |r cos[(co c + nw m )t + n0(t) + n |J 
n=0 


(62a) 
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(62b) 


s 2 (t) = exp^/3 sin|w m t + 0(tjj^ cos^co c t - /3 cosjo> m t + 0(tjjJ 

00 n 

= y (-l) n &T cos b + nw m )t + n0(t) + n |j 
n-0 ' L J 

S 3 (t) = exp j^/3 cos |w m t + 0(t)jj cos^ c t + /3 sin[w m t + 0(t)j| 

OO 

= ^ ^cos\(, c + n. m )t + ne(t)\ 
n=0 

S 4 (t) = exp |-/3 cos|w m t + 0(t)j| cos( Wc t - /3 sinjoJ m t + 0(t)jj 

00 n 

= ^ (-l) n ^-cosjjw c + nw m )t+ n0(t)] 


The power spectral density for any of the four signals, whether modulated or 
unmodulated, may be represented as 


S(o>) = ^ Sn(w) 


n=0 


where for the case of unmodulated subcarriers, 



and for the modulated case, 



where S mn (o>) is the spectral density of 

m n (t) = exp[]n0(t)j 


(62c) 


(6 2d) 


(63) 


(64) 


(65) 


( 66 ) 
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It is easily verified, using Parseval’s relation, that 


2tt 



S n (co) do> = 



( 67 ) 


for both equation (64) and equation (65). Thus, as is intuitively expected, the presence 
of angle modulation on the subcarrier does not change the total power in each spectral 
component, but just spreads the spectrum. 


Linear Product Demodulation 

The product demodulation of a sinusoidally phase-modulated SSB signal in the 
presence of noise is considered next. Figure 3 shows the model of the demodulator. 

Input Output 


s in (t) +n in (t) s 0 (t) + n 0 (t) 



s refW 


Figure 3.- Linear product demodulator. 


The input signal and noise terms are taken at the output of an ideal band -pass filter 
which drives an ideal product device. The input signal Si n (t) is assumed, ideally, to be 
undistorted and is taken as 

Sj n (t) = A exp (-/3 sin w m t)cos(cu c t + /3 cos «mt) (68) 

The noise is taken as ideally band-limited white Gaussian noise: 

n in(t) = x(t) cos w c t - y(t) sin u> c t (69) 

The reference signal s re f(t), which drives the other input port of the product device, is 
taken as 


s ref00 = -2 sin w c t 


(70) 
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Note that the reference signal is taken to be phase coherent with the input signal. This 
is not a loss of generality, since signals of the type dealt with here have a residual car- 
rier component of frequency u> c which may be detected and used to synchronize the* 
reference signal. Now the output signal s Q (t), the result of demodulating the first-order 
sideband, taken at the input to the band-pass output filter, is, neglecting terms of 
frequency 2w c , 


s Q (t) = A/3 cos co m t (71) 

Likewise, the output noise term n 0 (t), taken at the input to the band-pass output filter, 
is 


n Q (t) = y(t) (72) 

Now, the input signal s^ n (t) has only an upper sideband. Hence it is sensible to 
restrict the filter bandwidth Bi n to the upper sideband. As a matter of fact, the filter 
need not pass all the signal components. It need pass only the component at frequency 
oj c + o> m . In case a subcarrier system is considered, the filter should be wide enough to 
pass the modulated subcarrier without appreciable distortion. However, the filter need 
not pass the residual carrier (assuming it is detected in a separate channel) or higher 
order modulation components. 

Let nj n (t) be described by its spectral density N^ n (c<j). Since n in (t) is a real 
process, N in (co) is necessarily two-sided. That is Ni n (u;) exists for positive and neg- 
ative <x). Under the ideal filtering assumption, Ni n (co) has constant value, say ^ V* 
in the frequency neighborhood occupied by the signal component at frequency co c + co m . 
For all other frequencies, N in (co) is zero. Figure 4 gives Ni n (co) graphically. 


Nin(w) 

Mi,,! 








— 2 71 


I 
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9 

i 
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1 



-C0 c o w c w c + o) m 

Figure 4.- Noise spectral density Ni n (co). 


CO 


The input signal -to -noise ratio (SNR), the ratio of total modulated carrier power to 
total noise power, at the input to the product device is 
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(73) 


J m 


2 IqC 2 / 3 ) ^ a 2 

2 riBi 


Nin 1 „ OR- 
2 V “^m 


"in 


The output noise spectral density N 0 (co) is just the spectral density of y(t), or 
Ny(w). From appendix A it may be determined that N y (c«j) is equal to a constant, 

i r\, in the frequency domain of the signal component at frequency u> m and is zero 
elsewhere. Figure 5 shows N v (u>). 
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Figure 5.- Noise spectral density Ny(co). 

The output SNR at the output of the product device is computed for the noise 
existing in the output bandwidth B 0 . It is assumed that B 0 ^ Bj n . Then 


( Ml 


N ° h 2B o 


,2 3 2 
2 ??B 0 


.= = 


®in 

r ? 2 " 

A 2 B in 

e 2 1 

B 0 

Lio< 2 «J 

2 *? B in ~ B 0 

I 0 (2/3)_ 


3 in 


N. 


(74) 
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The function of /3 in brackets may be termed the "modulation gain" of this partic- 
ular modulation system. This function is the ratio of recoverable sideband power to total 
modulated carrier power given in equation (32). Thus the modulation gain factor has an 
absolute maximum of 0.475 for /3 = 1.29. Hence, it is seen that this type of modulation 
system is inherently a "low-index" system. 

In this section, analytic function theory has been applied to investigate the synthesis 
and detection in noise of single -sideband signals, phase -modulated by sinusoids. Mathe- 
matical expressions for recoverable signal power and for total modulated carrier power 
were derived for the case of modulation by a single sinusoid. It was shown that the power 
expressions are insensitive to angle modulation of the modulating sinusoid (subcarrier). 

It was also shown that for modulation by multiple sinusoids, the total modulated carrier 
power increases with increasing numbers of sinusoids, and that the amplitude of each is 
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independent of the others. It was determined that for linear product demodulation, this 
type of modulation system has highest efficiency for low modulation index (1.29). A 
synthesis procedure for signals of greater modulation efficiency was outlined. 

GAUSSIAN SINGLE -SIDEBAND FREQUENCY MODULATION 


In this section, SSB frequency modulation of a carrier by a Gaussian process is 
considered. Nonlinear (discriminator) detection of such a signal in the presence of 
noise is modeled mathematically for the high SNR case. Gaussian modulation is chosen 
since other message functions are often approximated by using a Gaussian process of 
suitable spectral shape. Frequency modulation is chosen as being representative of the 
"wide-band" angle -modulation process. 


Properties of the Modulated Carrier 

Since the signal treated here is frequency modulated, the message function is <£ s (t) 
where <£ s (t) is the phase function of the modulated carrier, as given in equation (10). 
Clearly, if $> s (t) is specified to be a Gaussian process, then 4> s (t) is also Gaussian, 
since the time derivative is linear. Thus, let 4> s (t) be defined as 

<Mt) = f (t) ( 75 > 


where f(t) is a Gaussian process, assumed weakly stationary and with zero mean. Let 
f(t) + jg(t) be an analytic Gaussian process. Then the analytic modulation function m(t) 
is given by equation (16) and repeated here for convenience: 


m(t) = exp 



ex p(j 0) 


where 0 is uniformly distributed on irj. Then the analytic signal is 


(76) 


^(t) = Am(t) exp(ju> c t) 


(77) 


and the autocorrelation function is 


R 1 //^ T ) = A. 2 ;t)| exp(jo> c r) (78) 

where A is the amplitude constant and <f> the characteristic function (ref. 9) given as 


^(l,j,-l,j;T) = E <exp 


jj[f(t+T) 


+ jg(t+T) - f(t) + 


jg(t)]j 


(79) 
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From equation (B22), R^(r) is developed as 


R^r) = A 2 exp ^Rg g (0) - %(0) + Rgg(T) + Rf f (r) + j [w c t + Rgf(r) - RfgfrjjJ 
= A 2 exp [2Rf f (r)] exp { j [“ c r + 2R g£ (r)] j 


( 80 ) 


It is interesting to note that the autocorrelation function for a double -sideband Gaussian 
angle -modulated carrier is obtainable from equation (80) as a special case. If g(t) = 0, 
the real part of the analytic signal of equation (77) becomes 


s(t) = A cos |u> c t + f(tjj 


(81) 


for which 



= exp (0)j exp[R ff (T)] cos u> c r 


(82) 


Now, for the SSB signal, taking the real parts of the various analytic expressions 

gives 



s(t) = A exp£-g(tj] cos£c<; c t + f(t) + 0] 

(83) 

and 

A 2 

R SS ( T ) = ~2~ 6XP [ 2R ff cos [-c T + 2R gf ( t 3 

(84) 

The "average power" of s(t), given by its variance a g 2 , is denoted P g and is 



2 

p s = P S 2 = R ss(°) = \ exp[2Rf f (og 

(85) 

since Rgf(r) 
function f(t). 

is an odd function. But Rff(0) is simply the variance o^ 2 of the phase 
Hence, 


P s = ^exp(2a f 2 ) 

(86) 


26 


Thus, equation (86) shows that the average signal power increases as the exponent of 
twice the message variance. This result is analogous to equation (30), which gives the 
average power for sinusoidal modulation. 

The Phase Process of Signal Plus Noise 

The frequency demodulator which is to be treated is sensitive not only to the modu- 
lated signal but also to the additive noise (Gaussian) in the signal channel. The demodula 
tor detects the time derivative of what may be called the Equivalent phase TT of the signal 
plus channel noise. Appendix C treats the equivalent -phase concept in detail. 

Let the noise n(t) and signal s(t) be given in the standard forms of equations (4) 
and (10), respectively. It is possible, in general, to represent the sum r(t) of the 
frequency -modulated signal s(t) plus noise n(t) as 

r (t) = s(t) + n(t) = R(t) cos £u> c t + <p( t) + 6>(t) + (87) 

where R(t) is nonnegative, t) is the signal phase function, <9(t) is the perturbation 
of the equivalent phase due to additive noise, and 6 is a uniformly distributed constant. 

In the case at hand, 


A(t) = A exp [-g(t)] 
0(t) = f(t) 


( 88 ) 


where f(t) and g(t) are jointly Gaussian and f(t) + jg(t) is analytic. By Substituting 
in equation (Cll), the low-noise approximation to 0(t) is obtained as 


0(t) ~ ^ exp[g(t)] jy(t) cos [*3 - x(t) sin[f(tj]l 
Noting equation (B23), R qq(t) is obtained from equation (C18) as 

R 0$( T ) = | Re {^~2[ R zz( T ) J 

where the characteristic function is evaluated in equation (B22). Then 

r 00( t ) = ex P pggC 3 ) " Rff(°) + R gg( T ) + Rff ( T j] 

^ R xx(^) cos|2Rgf ( t)J - Ryjj(T) sin [ 2 %w]} 


A" 

x 


(89) 


(90) 


(91) 
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Equation (91) is not the most compact form but is given so that a specialization may be 
made to classical double -sideband FM. Notice that the autocorrelation function for the 
phase-noise process of a classical FM signal in noise is obtained by letting g(t) = 0 as 
in equation (81). Then, from equation (91), 

r 06 >( t ) = 4 > ex P [-Rff(O)] Rxx( T ) ex P [%( T )] ( 92 ) 

Now, it is common practice in classical FM to assume that, in the low-noise case, 
the phase noise is simply 


m = 


(93) 


which gives 


R ee ( T ) - a 2 R yy^ ( 94 ^ 

This result is not rigorously correct. However, for high-index modulation, exp Rff(r) 
is much broader than R xx (r). Hence an approximate result, good for high-index modu- 
lation, is 


R 06 )(t) * ^ exp [-%(°3 RxxW ex P[%(°)] 

A 

= ~^2 -^xx( r ) 

A 2 

= ~2 %y( T ^ (95) 

For SSB-FM, equation (91) may be written more compactly as 

H 00 (r) = ex P [ 2R ff( T )] (Rxx^ cos [2RffW] " Rxx(t) sin [ 2 % (r)] j (96) 

As in classical FM, if the signal channel bandwidth is sufficiently wide compared with the 
bandwidth of the modulation function f(t), then equation (96) reduces to 

R ^( T > *72 exp ^ R ff(°)] RxxW (97) 

A. 
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Limiter -Di sc riminator Demodulati on 

Now the response of an idealized limiter -discriminator to an SSB-FM signal plus 
band-limited white Gaussian noise will be determined. The model of the demodulator is 
given in figure 6, where all filters are taken to be ideal. 

Input s inW +n in(t) 



s d(t) + nd(t) s 0 (t) + n 0 (t) Output 

Figure 6.- Demodulator model. 

Since the signal Sj n (t) at the input to the band-pass filter has only an upper side- 
band, it is sensible to use a filter which passes only the upper sideband. Then the noise 
process nj n (t) at the input to the ideal limiter also has only an upper sideband, refer- 
enced to oj c , the signal carrier frequency. The band-pass filter at the output of the 
limiter is zonal. That is, the filter passes only the first zone output of the limiter, 
which is a sinusoidal carrier of frequency u> c , having only phase modulation (ref. 15). 

Then the signal A(t) is represented as 

A(t) = cos £u> c t + f(t) + 0(t) + (98) 

where f(t) is Gaussian and 0(t) is the phase-noise perturbation. 

Since the objective is to determine the best performance of the limiter -discriminator 
in demodulating SSB-FM in the presence of noise, the low-noise approximation for 9( t) 
is made. Also, it is assumed that the input bandwidth B| n is sufficiently large with 
respect to the spectral width of f(t) so that Rqq(t) is well approximated by 
equation (97). 

The signal and noise outputs of the discriminator, s^t) and n^(t), respectively, 
are given by 
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( 99 ) 


s d(t) = ~5(t)] = f(t) 


n d (t) = 


The output filter is of the low -pass type, under the assumption that the message process 
f(t) is low-pass. The bandwidth of the filter B 0 is chosen so as to pass f(t), or most 
of the power in the spectrum of f(t), while not passing any more noise than necessary. 
Thus the output signal power S Q is essentially 


S o = Rs d s d (0)=-^[%W]|^ 0 ^-Ri;(0) 


( 100 ) 


To determine the output noise power N 0 requires first the determination of the noise 
spectral density, and hence the autocorrelation function of n^(t): 


^d 1 ^^ “ ~ ^2 [^•06 | ( t )J " ~ R 0d( T ) ~ ” ~2 exp [ 2R ff(°)] r xx( t ) 
Thus, the spectral density N n( ^(w) of n^(t) is 


( 101 ) 


N nd M = exp[2R ff (0)] co 2 N x (co) 


( 102 ) 


where N x (w) is the spectral density of x(t). If the input spectral density N in (o>) is 

equal to a constant, say ~ r\ over the bandwidth Bi n , then from appendix A it may be 

seen that N x (w) is also constant over the bandwidth B Q and is also equal to ^ rj. 

z 


Then 




N n d ^) = exp [ 2R ff(°)] 


(103) 


The output noise N 0 in the bandwidth B 0 is just 


N 0 = — f 


27tB 0 12^ r -.9 / \ * 

N na Hd» = £ fj- exp [2%(0>] !(2,B 0 ) 


2?r J-27TB r 


A2 


(104) 


The output SNR is then 
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(105) 


So_ = ~ R ff^ 

N ° 

To evaluate equation (105) and relate S 0 /N 0 to the input SNR, a spectral density 
for f(t) is chosen. Let the spectral density for f(t) be ideal low-pass as shown in 


Sff(co) 



(jO 


-Cdf 0 Wf 


S|f(co) 



figure 7. Then f(t) is also low-pass. The autocorrelation function is then 

Qj sin (Ujt 


Rff(r) = 


7T T 


where Q f is the amplitude of the spectral density and 

-%(0) = = -f- Rff(0) 

Obviously, the optimum choice for output bandwidth in this idealized case is 


Wj = 27 tB 0 


(106) 


(107) 


(108) 


Then, for an ideal flat-output filter. 


31 



S Q 

No i r] 

%g> 2 ^ [ 2R ff (0) ] 

2 w f exp[4R ff (0)] 2 — — 770 J in 

2ir 2 in 


But by defining 


w in = 


( 110 ) 


it is seen that 






( 111 ) 


where S^ n /N^ n is the input SNR in bandwidth Bi n . Hence 

|=2^ % (0)exp[-4 % (0j]^. (112, 

Since Rff(0) is the mean-squared value, or variance, of the phase-function f(t), equa- 
tion (112) shows that for high-index modulation, SSB-FM delivers modulation loss rather 
than modulation gain. 


CONCLUSIONS 

This analysis produced two types of results. First, it showed that mathematical 
investigation of complicated single- sideband angle -modulation and demodulation problems 
is made tractable through use of the theory of analytic functions. Second, it gave tech- 
nical results concerning the communications efficiency of two SSB angle -modulated sig- 
nals. For the case of high-index SSB frequency modulation by a Gaussian random pro- 
cess, the results were negative. Investigation showed that even in the low -noise case, 
limiter -discriminator detection of such a signal produces modulation loss, rather than 
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gain. On the other hand, it was determined that low-index SSB phase-modulation, using 
subcarriers and linear product demodulation, is theoretically feasible. In fact, a synthe- 
sis procedure for such a signal with high modulation efficiency was outlined. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., May 14, 1969. 
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APPENDIX A 


SPECTRAL RELATIONS FOR BAND -PASS NOISE AND SIGNALS 


Noise Processes 

In computational applications, the spectral density of a band-pass process n(t) is 
usually given. When n(t) is expressed as 


n(t) = x(t) cos w c t - y(t) sin co c t 


(Al) 


the spectral densities of x(t) and y(t) must be determined. This determination may 
be made in general for a noise process of arbitrary sideband structure as follows. 

Let n(t) be the real part of an analytic process ^(t), where 


v(t) = z(t) exp(jco c t) 
z(t) = x(t) + jy(t) 


(A2) 


It is assumed without loss of generality that n(t), x(t), y(t), z(t), and v(t) are zero- 

mean weakly stationary processes. The autocorrelation fimction of z(t) is 


R-ZZ (^") ~ ^3Qf(^") Ryy( T ) ■*" 1 ^Ryx( T ) “ Rxy("^ 

Now, it is necessary for the weak stationarity of n(t) that 

Rxx( T ) = %y( T ) 

Ryx( T ) = -Byx( -T ) 


(A3) 


(A4) 


Thus, 

Rzz( T ) = 2 [Hxx( T ) + i-^yx( T )j 

R w (t) = R Z z( T ) exp(jw c rj > (A5) 

= 2 |Rnn( T ) + jRnn( T )] 

where the caret denotes Hilbert transform. The last equality of equations (A5) follows 
from the analyticity of u(t). 
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APPENDIX A - Continued 

The spectral densities are defined as the Fourier transforms of the respective 
autocorrelation functions and are given as 


S zz (u>) = 2 [S^co) + jS yx (co)] 


S w (w+u c ) 


S vv (u) = 2 [l + sgn(w)] S nn (cn) 


The second equality of equations (A6) follows from multiplying Rj,j,(t) by exp jjju> c Tj 
and applying the convolution theorem. 

Now by a result of Fourier transform theory (ref. 16), 


— ^ ” jSyx(uO" 


By addition and subtraction of equations (A5) and (A7), it follows that 


SxxM = 


®zz( w ) ®zz(' w ) 


S VX M = i 


S zz (-w) - S zz (u>) 


From equations (A6) 


S zz (w) = 2 |l + sgn(u>+w c )J S nn (u>+w c ) 
Szz(" w ) = 2 [l + sgn(-w+w c )j S nn (-cu+w c ) 


Thus 


®xx( w ) - ®yy( w ) 


2 [l + sgn(oH-u> c )] S nn (uH-o> c ) + 2[l + sgn (-UH-o> c j] Sn n (-w+w c ) 


[l + sgn(oH-u> c j] S nn ( om-c<; c ) + [jL - sgn(w-u> c )] Snn( Ct,-w c) 


(A10) 
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APPENDIX A - Continued 
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-2co c -co c 0 co c 2co c 

Figure Al.- Spectral relations for noise. 

by the "evenness" of S nn (co) and the "oddness" of sgn(co). Equation (A10) can be 
written in terms of the unit step function U(co) as 

Sxx(co) = Syy(co) = S nn (om-u> c ) U (cu c +o>) + S nn (co-co c )u(co c -co) (All) 

An interesting observation about the correlation of x(t) with y(t) may be made 
by consideration of the cross spectral density Syx(co). From equations (A8) and (A9), 

SyxM = j|s n n( w - w c) u ( w c- w ) ~ s nn( w + w c) u ( a; c +w )] ( A12 ) 

Since S nn (co) is real, S yx (co) is purely imaginary, a consequence of the "oddness" 
of Ryx(r). Now, although it may not be obvious from equation (A12), if S nn (co) is 
locally symmetric about o> c , then S yx M is identically zero. Thus Ry X (T) must be 
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identically zero. Hence x(t) and y(t) are orthogonal. Since x(t) and y(t) are 
assumed to be zero-mean they are also uncorrelated. To illustrate equations (All) 
and (A12) for a process which is locally symmetric about co c , figure A1 is given. 

Signal Processes 

For signal processes, the computational problem is generally the inverse of the 
problem encountered with noise. That is, the spectra associated with the modulation 
function are generally known and the spectrum of the modulated carrier is to be derived. 

Let a real band-pass process s(t) be given as the real part of an analytic 
process i//(t) by 


s(t) = A(t) cos [c c t + <p{ tj] 

^(t) = m(t) exp(jw c t) 

The "modulation process" m(t) may be represented, as for noise processes, by 

m (t) = x(t) + jy(t) 


(A13) 


x(t) = A(t) cos <p(t) 


y(t) = A(t) sin <p(t) 




(A14) 


As with noise processes, it is required that s(t), x(t), y(t), m(t), and \p(t) be zero- 

mean, weakly stationary processes. 

Because of the weak stationarity of s(t), equations (A4) hold for the present case 
also. Then 


Also, 


Thus 


Rmm( T ) = 2 [Rxx( T ) + jRyx( T )] 

R^( T ) = 2[R SS ( T ) + jR S s( T )] = R mm( T ) exp(jo) c r) 

S^M = 2(1 + sgn cu) S ss (cu) = S mm (w-o> c ) 
S^(-co) = 2(1 - sgn co) S ss (cu) = S mm (co+o> c ) 


(A15) 


(A 16) 


(A17) 
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Sxx( W ) - Syy (Oj) 


0 

jSyxM 

y\ 


Ki 

^XX^+^c) + Sxx( C0-a, c) 




1 1 

T 

-to c 

0 

J jSyxC^-^c) " Sy X (^+^c)] 

co c 


0 

Sss(w) 

l^c 

-o) c 

0 

w c 


CO 


CO 


CO 


CO 


Figure A2.- Spectral relations for signal. 


It follows that 


Then 


2(1 + sgn w) S ss (o>) = 2js xx (w-a; c ) + jS yx (w-a> c Jj 
2(1 - sgn w) S ss (u>) = 2[s xx (w+w c ) - jSy X (a>+a> c )] 


2S ss M = Sxx( w - w c) + j s yx( w " w c); 

2S ss (o)) = - jSy X ^a>+0) c 'j ; 


(a) > 0) 

(w < 0) 


Because of the analyticity of 


(A18) 


(A 19) 


s mn» = °; (M = w c) 


(A20) 
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Thus 


’ss^) - 2 ^ ^xx( w w c) + ®xx( w+Ct) cjj + j|®yx( c<J w c) “ | (A21) 


It should be noted that because Ry X (r) is real and odd, S yx (co) is odd and imag- 
inary. Hence S ss (co) is real and even, as required of spectral densities of real pro- 
cesses. To illustrate equation (A21), figure A2 is given. 

In computational work it may be easier to compute S ss (cu) directly from equa- 
tions (A17) as 


^ss^) " + ^mm( co+co c) 


(A22) 
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STATIONARITY CONDITIONS FOR STOCHASTIC 
SINGLE -SIDEBAND ANGLE MODULATION 

The General Case 

The analytic form of an SSB angle -modulated signal s(t) may be written as 

i//(t) = m(t) exp(jw c t) 

For s(t), the real part of i//(t), to be weakly stationary requires that 

E {m(t)} = E ^m (t+r) m (t)} = 0 


(Bl) 


(B2) 


Now consider m(t) of the form 

m(t) = exp{ j 


j[f(t) + jg(t)]j 


(B3) 


where f(t) + jg(t) is analytic. Because the exponential function is entire it follows that 
m(t) is analytic, and hence i^(t) represents a signal having only an upper sideband. 
The conditions of equation (B2) are then written 

E {m(t)} s E<jexp ^j[f(t) + jg(t)]|^j> 4 $(l,j) 4 0 (B4) 

E {m(t+T) m(t)} s E^exp|^j[f(t+T) + jg(t+T) + f(t) + 4 $(l,j,l,j;r) 4 0 (B5) 

where <f> denotes characteristic function. 

Thus, the stationarity requirement of equation (B2) is written as 


*(i*3) = = o 


(B6) 


where the ordering of the stochastic variables is understood to be that implicit in equa- 
tions (B4) and (B5). 


The Gaussian Case 


Now let f(t) + jg(t) be a complex Gaussian process, 
individually and jointly Gaussian. Assuming that f(t) and 


Then f(t) and g(t) are 
g(t) are zero-mean, they 
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are also individually and jointly stationary in the wide sense. The characteristic func- 
tions may be expressed (ref. 9) as 

= E fexp^WjZi + o> 2 z 2 + ... + w n z n 



s exp 


n n 


sll 

j=l k=l 


(B7) 


where 


"jk 4 E{z jZk } 


The terms 


[tv*] = 


re easily tabulated 
~R ff (0) jRgf(O) 

jRgf (0) -R gg (0) 

R ff (T) jRgf(r) 

jRj g (T) -R gg (r) 


in matrix form as 


R ff (r) 

jRfg(T) 

jR g f( T ) 

-V T > 

R ff (0) 

jRgf(o) 

jRgf(0) 

- R gg<°> 


for the characteristic function <3>(l,j,l,j;r). Also, 



3R fg (0) 

-Rgg(o) 


(B8) 


(B9) 


(BIO) 


for the characteristic function $(l,j). 

Now, from equations (B7) and (B9), 


*(l,j,l,j;r) = exp|Rgg(0) - Rf f (0) + R gg (r) - R ff (r) - j[2Rgf(0) + Rf g (r) + Rgf(r)]J (BU) 


From equations (B7) and (BIO) 


*(l,j) = exp^i|Rgg(0) - R ff (0)] - jRgf(0)j 


(B12) 
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It is obvious from equations (Bll) and (B12) that the stationarity conditions of equa- 
tion (B6) are not satisfied for m(t) of the form of equation (B3), where f(t) + jg(t) is 
Gaussian. 

Consider a modified modulation function m'(t) such that 

m'(t) = m(t) exp(j 0) (B13) 

where 0 is a random constant uniformly distributed in Q-tt, tt] and statistically inde- 
pendent of m(t). Then 

E{m '(t)} = E{m(t) exp(j0)} = E{m(t)} E{exp(j0)} (B14) 

Now 

E{m(t)} = <3?(l,j) ^0 (B15) 

which implies 

E{m ’(t)} = 0 if and only if E{exp(j0)} = 0 (B16) 

Likewise 

E{m'(t+r)m'(t)} = E{exp(j2 0)} (B17) 

which implies 

E{m'(t+T) m ' (t)} = 0 if and only if E{exp(j20)} = 0 (B18) 

It is readily verified that for 0 uniformly distributed in {-tt, 7 7 ], 

E {exp (j 0)} = E (exp (j 2 0)} = 0 (B 19) 


Hence, 


m’(t) = exp 



exp(j 0) 


(B20) 


satisfies equation (B2), the stationarity requirement on s(t), where f(t) + jg(t) is ana- 
lytic and Gaussian. 
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Expansion of the characteristic function 4?(l,j,-l,j;r) is required for the 


["jVjkl - 


WjCOkPjk 

terms for this characteristic fun( 

Rff (0) 

jRfg(O) 

-Rff(r) 

jRfg(T) 

jRgf(O) 

-Rgg(°) 

-jRgf(-r) 

- R gg (T) 

-R ff (r) 

-jRgf(r) 

Rf f (0) 

-jRfgfO) 

iRfgW 

- R gg< T > 

-jRt g (0) . 

- R gg (0) 


It follows that 


= exp |ltgg(0) - ^(0) + R gg (T) + Rf f (T) + j[Rgf(r) - Rfg(T)jJ 
= exp ^2 RffCr) + jR gf (T)~ J 


Also, it is easily shown that 




text. 


(B21) 


(B22) 


(B23) 
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THE EQUIVALENT PHASE PROCESS OF THE SUM OF AN 
ANGLE-MODULATED CARRIER AND GAUSSIAN NOISE 


Let an angle -modulated sinusoid be represented in standard form as 

s(t) = A(t) cos[o> c t + $(t)J 

where <fi( t) represents the angle modulation and A(t) is a nonnegative "envelope 
function.” Let an additive noise process be given in standard form as 

n(t) = x(t) cos 0 ) c t - y(t) sin cu c t 

The analytic forms xf/( t) and u(t), of s(t) and n(t), respectively, are 

i//(t) = A(t) exp|j0(t)J exp(jo; c t) 

Kt) = z(t) exp(jw c t); |z(t) = x(t) + jy(t)J 


(Cl) 


(C 2) 


(C3) 


Now, it is desirable to represent the sum of s(t) and n(t) as a single modulated 
sinusoid r(t) of the form 


r(t) = s(t) + n(t) = R(t) cosjoL> c t + <p( t) + 0(t)j 


(C4) 


in the manner of Rice (ref. 17). In this representation, R(t) is a nonnegative envelope 
function, cf>( t) is the angle modulation of the original signal, and 9( t) is a perturbation 
due to the additive noise n(t). The analytic form p(t) of r(t) is 


p(t) = R(t) exp|j0(t)j exp[j0(tj] exp^jw c tj 
Now, the sum of ij/( t) and v(t) is 

H t) + 2 /(t) = ^A(t) expjj0(t)| + z(t)j exp(jw c t) 

Hence, equating equations (C5) and (C6) and dividing by exp(jo> c t^ gives 

R(t) expjj0(t)J exp|j 0(t)j = A(t) exp|j0(t)j + z(t) 


(C5) 


(C6) 


(C7) 
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and 


R(t) exp[j<9(t)] = A(t) + z(t) exp |^j 0(t>] 
Equating real and imaginary parts of equation (C8) gives 


(C8) 


R(t) cos 0(t) = A(t) + x(t) cos <p( t) + y(t) sin (p(t) 
R(t) sin 0(t) = y(t) cos <p( t) - x(t) sin cp(t) 


(C9) 


It follows that R(t) is nonnegative and that 


0(t) = tan" 


y(t) cos 4>(t) - x(t) sin <p(t) 


A(t) + x(t) cos <p( t) + y(t) sin <$>{ t) 


(CIO) 


It is seen that, in general, 0(t), the perturbation of the sum signal r(t), is a non- 
simple function of the noise and modulation components. However, an approximation to 
0(t) may be obtained for the case where the mean-squared values of x(t) and y(t) are 
much less than the mean-squared value of A(t). The usual "small angle" approximation 
is 


6>(t) ~ y(t) COS ft® . - x(t) Sin ft® ; |A(t) ± 61 (Cll) 

A(t) A(t) L J 

It is assumed that in nontrivial cases the probability of A(t) = 0 is zero, and hence 
0(t) is given for small x(t) and y(t) by equation (Cll) almost all the time. 

Having now an approximate model for the "phase-noise" component of r(t), it is 
next desirable to develop the autocorrelation function of 0(t) with a view toward 
obtaining the spectral density: 

0( t+T ) 0 (t) « | y(t+r) cos 0(t+T) _ x(t+r) sin 0(t+T)~ | f y(t) cos 4>(t) _ x(t) sin <j)(t)~ | 

L A(t+r) A(t+r) JL A(t) A(t) J 

_ f y(t+T) y(t) cos 0(t+r) cos <p(t) t x(t+r) x(t) sin </>(t+r) sin 0(t) 

L A(t+r) A(t) A(t+r) A(t) 


x(t+r) y(t) sin <p(t+r) cos 4>(t) _ y(t+T) x(t) cos 0(t+r) sin <p(t) 
A(t+r) A(t) A(t+r) A(t) 


(C12) 


45 



APPENDIX C - Concluded 


Now, assuming that x(t) and y(t) are both statistically independent of A(t) and 
4>(t), R ee (r) is 

R ed (r) = E (e(t+r) 0(t)} 

= R vv (t) E (c° s ^(t + T) cos 0(t)] + r (t) e fm M + j) sin 0(tfl 
Ty l A(t+r) A(t) J ** I A(t+r) A(t) / 


- Rx V (t) cos 0(tA _ p (t) e ( cos <p (t+ t) sin (cl3) 

y I A(t + r) A(t) J yX l A(t+r) A(t) J 


Because n(t) is Gaussian, 

R 06 ,( t ) 


_ p / T \ E fcos 0(t+r) cos 0(t) + sin 4>(t+ r) sin <p(t)\ 
^ l A(t+r) A(t) J 


+ R (t) e f sin < ft( t+T ) cos 0(1) ~ cos sin 0(t)"l 

L A(t+r) A(t) J 

= (T) E f c°?&(t) - »(t+rjn _ E r S in[0(t) - 0(t» T a~i 

l A(t+r) A(t) J l A(t+r) A(t) j 


Now, since A(t) is nonnegative, it may be represented as 

A(t) ^ A expj^-o'(t)J 

Then, denoting 

^ E /exp|j[-0(t+T) - ja(t+T) + (p(t) - ja-(t)jj 


it follows that since 


then 


Rzz( T ) — ^ [Rxx( t ) + jRyx( T )] 


R 00 (t) = | Re (^2 r zz(t) 


(C15) 


(Cl 6) 


(C17) 


(C18) 
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